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The external electromagnetic field changes the 
states of system. The zeroth, first and second or- 
der corrections in wave function required for con- 
ductivity and Hall mobility are listed. In general, 
there are several types of processes contributing to 
a given transport coefficient. In each type of trans- 
port process, for a group of given electronic states, 
the electronic and vibrational contributions are fac- 
torized. 



The vibrational degrees of freedom can be in- 
tegrated out, and result in a time integral. The 
time integrals depend on the electronic state, tem- 
perature and the frequency of external field. The 
time integrals for conductivity and Hall mobility 
are given. The time integrals can be estimated by 
asymptotic analysis. We sketch how to implement 
the present results in ab initio code. 



1 State and the changes in state caused by exter- 
nal field To compute the conductivity, one needs t^(°) 
and tp^\ while to calculate Hall mobility one needs 
ip(°\ and xp( 2 \ cf.[text, Eq.(17)]. As discussed in 
Sec. 4, the lowest order self -consistent approximation 
requires computing all these quantities to order S° 
andS 1 . 

1.1 Initial state is a localized state By the proce- 
dure sketched in the end of Sec. 4, to first order in 
small parameters / and ]' the state of system is 
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To first order in external field, the change in state 
i/>W to S° and S is 
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Because the carriers obey Fermi statistics, in Eqs. (|ll2[> , 
a factor 1 — f(E Al ) must be included before compo- 
nent (pA 2 of the final state, factor 1 — /(EgJ must be 
inserted before component £b 1 of the final state. 

The full expression for xp^ is too long to repro- 
duce. Supposing there exist only localized states, to 
order J° and J 1 , tp^ is: 



e it /,f {S^ } -s[ Na} )/h 
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1.2 Initial state is extended state If the initial 
state is \B • • • N a • • • ), the state of system at time t 
without the external field is 

^(°)(o = ?B4 N " } e- i <" }/R 
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The last term is an order K 2 contribution, and is use- 
ful only to calculate conductivity from EE transition. 
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To first order in the external field, the change in state 
caused by the external field is 
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it'(£^-£^)/h 



U"(£if-£^)/n_ 

In Eqs. (|4l5|) , a factor [1 — f{E Al )] must inserted before 
component (p Al of the final state, a factor [1 — /(Eg )] 
must be included before component £g of the final 
state. 

2 Sum over final phonon states and average over 
initial phonon state For processes involving only / 
and ]' , one first calculates inner products to second 
order in the static displacements (Q A , Q Al , 9 Al , A3 ). 
Then one multiplies the inner products, neglects the 
terms of third and higher order in (0 A , 6 Al , Al , A3 ) 
and the terms which involve more than one phonon 
changed in one mode. One can first compute the sum 
over phonon final phonon states and average over 
initial phonon for a single mode, and then multiply 
all modes together. 

For processes involving e-ph interaction K' or 
K, the matrix element of K' or K has shape of 

[Loc> /n", n' ] [ Pn", n'}- The sum over final phonon 

oc" oc' oc' 

states and average initial phonon state takes the form: 

E I1W E IYIsn»,nJ e 



■Not- * 



,. N »... oc 



E/n"N' EI PN"M f JlTI h K,Nj' ( 6 ) 

oc' oc(^oc') ol 

Using the induction method, one can see that Eq.© 
is equivalent to: 

EEiS^N", fN'^N' t gN"„N ttt hN' f ,N, . 
oc oc oc oc oc ft oc' % 1 

a r En',N", Pn", ,N' f &N",,N i ^N r f ,N / 

u lv oc' oc' oc' oc' & oc" & 



IK E PK',N'„gN'J,N x h K,Nj}- 



(7) 



In the first square bracket in Eq.((7j, each term in the 
denominator is a pure number (does not include Kro- 
necker delta signs). One can expand each fraction into 
a power series in the static displacements (0 A , 9 Al , 
9 Al , 6 A3 ), neglecting terms of third and higher order 
in (6 A , 6 Al , Q Al , A3 ) and the terms in which more 
than one phonon changed in one mode. The remain- 
ing calculation is straightforward. 

Taking Fig.5(a) as example, according to the rule 
in tables 2 and 3, the contribution to current density 
is 

(-N e he/m) £ f(E A )YlP(N 0i ) 

A — Noc— oc 



E [i-/(£a 2 )] E 

a 2 -K- b v -K- 
Q" 1 / dr(^V r ^ 2 -^ 2 V r ^)(Yif>|Yf^) 

e it(£if-£^)/n 

^4f } i< 4 )7_l^'^£ w (^)]* (8) 

Expanding Y^"^ to second order in the origin shift 
9 A , and using the orthogonality of harmonic oscilla- 
tor wave functions, one has 

<4f rt N ' } >=n*NLN.< < 9 > 

oc 

where 

+ (tf)l(f ) 1/2 ^- lX - (^) 1/2 <W X ]- 



Similarly, one has 

oc 

where 
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The matrix element of the e-ph interaction isllTI 

Mf } i^ 2B i4f > = E/ W i n pn» 
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+^ 2 [(f ) 1/2 ^-iX'"( 



N a + l^l/ 2A i 



A simple and lengthy calculation shows that 



ESn"„n,/ lsN "'^' E/n"„n', V,N,/ 
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E^n» N/" isN ^ E Pn-X^n/"^ (17) 



A\2 



, N c + 1 
H ^ — e 
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Combine fL6)) and fL7|) , we have 

LN , f N"jN" f X,8N" / ,N, h N , / ,N, e 

DC DC DC DC DC & DC & 



v-> . . _ r, -is'N' ,cv a /-isN",cv„ f 



■N' , N", VN", ,N' , SN". ,N„, n N r „N„r 6 
a' a' a' a' a' a a' a 



q^2 



(18) 

Making use of Eqs.0, ([17)) and ([18}, one reaches 
Eq.©. 

To fulfil integrating out the vibrational degrees of 
freedom, two approximations are involved [2]. When 

computing inner product (Y^^Y^^) and matrix 

element (Y^ } \K' AiB |sjf" } ) etc , we neglected (1) the 

terms which are (0^) 3 and higher order; and (2) the 
terms in which two or more phonons are changed in 
one mode. When temperature is close to the melting 
point, the present results may be only qualitative. 



3 Time integrals for conductivity 

3.1 Localized state as initial state Fig.2(a): 

Ia 1 a±(w) = f° dse ±iws e~ i<E ' A i~ E ' A)/h (19) 

J — 00 

exp{ — 2^ [coth —— — (1 — cos cOccS) — i sin cv^s] } 

a 

Fig.2(b): 

Ib,a±= [ dse ±iws e- is{EB i- EA],h (20) 

J — 00 

exp(-E{^oth^(^) 2 

a 

— -(0^ ) [coth — cossojoc — zsmso^Jj), 
Fig.3(a): 

Ia 3 a iA ± = [° ds f° ds'e^' (21) 

7—00 J — oo 



e i S (E Al -E A )/H-i S '(E A3 -E A )/h gxp( _ £ { 1 CQth PhWa 

DC ^ 2 

4 ^ - « ) - e « 1 ) ( c ° th ^ c ° s ^ - * sin ) 



[coth ffi^ cos(s — s')^ + z sin(s — s')a; a ] 

+ 5 (0a 1 ~ Of) (^ 3 - ^ ) (coth ^ cos scv. + z sin so;,) }) 
Fig.3(b): 



h 1 A 3 A± = ds ds'e 

J —oo J — 00 



r n ±iws' 



(22) 



e is(E A3 -E A )/h-is'(E h -E A )/h 

exp(- E(l ^th ^[(^)2 + ( ^ )2 + _ ^ )2] 

a 

— 2 a a ( co 2 — C ° SS ^ _ zsms 
+ ^0* 3 (0* 3 - 0« ) (coth cos scooc + z sin sa; a ) 

(^ 3 - 0^) [coth ^ cos(s - sf)a> a + i sin(s - s f )cv a )}) 
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Fig.3(c): 



a 2 a± = [° ds f° ds'e ±iws ' (23) 

J —CO J — CO 



I 



e is(E B3 -E A )/h-is'(E A2 -E A )/h 

ex P (- U\ coth ^ [{& - + + (^) 2 ] 

+ ^* 2 (^ 2 - 0« ) (coth cos s'w* - i sins' Wo,) 



"2^ot\yoL -Vol) 



(coth p - - cos(s — s')^ + i sin(s — s')a; a ) 

— ^0^0^ 2 (coth ffi^ cossa^ + zsinsa^)}) 
Fig.3(d): 

I B3BlA± = /° /° ds'e W (24) 

J — CO 7—00 



[coth — - — cos sa^ — z sin sco^ }. 

Fig.4(b): 

Q 2 ^ 1 A±(^T)= /° d se -»( E ^- E Ai)/» (27) 

J— 00 

/° ds'^^e- i ^ +8 ^ EA i- E ^ /h W A2AlA (s,8' r T). 

J — CO 

Fig.4(c): 

Qib 2 W^T) = /° dse ±isw e- is ^2- E ^ /h (28) 

J —00 

J —CO 

where 

W B2 a 1 a(s,s , ,F)= (29) 
exp{- \ Ecoth ^[(ctf) 2 + (^) 2 + (tf * - O 2 ] 

^ DC Z 

+ J E « 1 ( « 1 - e « ) I coth ^ cos - i sin s '"«] 



[coth — coth — - — cos(s — s )w a — i sin(s — s Ja> a J}. 

Fig.4(a): 

Q 1A2AlA± (u,T) = f ds[° ds ' e ^ e -HE A2 -E Al )/H 

J —CO J —CO 

(25) 

e-^' +s ^- E ^ /h W A2AlA (s,s',T), 

where 

W A2AlA (s,s',T) = eX p{-iEcoth^ (26) 

[(0« 2 - + (^ 2 - e* 1 ) 2 + (tf 1 - et) 2 } 



-\ E(^ ! " #)#[«>th ^ cos( S + s')a; a - zsin(s + s> a ] 

a 

+ ^ I] a a 1 [ Coth ^—y^ COS Sa; ^ - i Sin SCVcc] } . 

Fig.4(d): 

Q2B 2 A lA ±(co,T) = f° dse- is ^- E ^/ n (30) 

J —CO 

f° ds ' e ±i(s+s')u, e -i{s+s')(E Al -E A )/n WB ^ SrS ^ T y 

J —CO 

Fig.5(a): 



Q 



1 



[coth i — ^ — cos s^-i sin s co a \ 



K' 

1A 2 B 1 A± 




[a>,T)= f° dse~ is{EA 2- E ^ /h (31) 

J — oo 



/U 
-oo 



.±t(s+s')w -i(s+s')(E Bl -E A )/ft 



+jE(^-^)(^ 2 -^) 



[coth p - - cos(s + s')cOpc — i sin(s + s')a; a ] 



V^aW 2 

{0$ (i coth sin s'ay - cos s'av ) - 

+ ( e ft - «> )(« coth sin sa V - cos sa V ) }] 

exp{-i E^th ^[(^ - ^) 2 + (^) 2 + (^) 2 ] 
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— 2^ ~% — [coth ^ — cos s^-isins o^J 

OL Z 

[coth cos(s + s')cOcc — i sin(s + s')^] 



Fig.5(b): 

Q2A 2Bl A±Kr)= /° d Sg ±&w e - & ( E ^- E Bi) /s (32) 

J — 00 



/° ds / e- i ( s+s ')( EB i-^) / *W A2 B 1 A(s,s' / T) / 

J — CO 



where 



(33) 



ex P {- \ Ecoth ^[(^ - ^) 2 + (^) 2 + (^) 2 ] 

^ OL Z 

+ - 2^ Vol ®ol [COth ' — - — COS S €0% — I Sin S CO^ 



[coth p - - cos(s + s')^ — i sin(s + s')a; a ] 
+ ^ E(^ 2 - ^)^ 2 [coth ^ cosset - i sin so; J}. 

Z OL 

Fig.5(c): 

Qib 2Bi a±(^T)= /° d se -»( E %- E B 1 )/» (34) 

.7 —00 

J— 00 



[(coth - — - — cos s av — z sin s ay j 
+ (z coth — - — sin say — cos say ) \ 



exp{-i Ecoth ^(O 2 



+^E(^) 2 



[coth ^ g cos(s + s')a> a — /sin(s + s')co a ]}. 



Fig.5(d): 

Q2B 2Bl A±(co r T) = f° dse ±isu e- is ^2- E ^ /n (35) 

.7 —oo 

/° ds'e-^'^i-^expi- 1 - Ecoth ^(0 a A ) 2 

J— 00 2 ^ 2 

+ ^ E ) 2 [ Coth COS ( S + S ^ ~ 1 sin ( S + S ^]}' 

3.2 Extended state as initial state Fig.6(a): 



j — ( 



Ia 2 b± = I dse ±iws e- is{EA 2- E ^ ,h (36) 



exp{-icoth^(^) 2 

+ 2^a; ) [COth — COS SCOqc — zsinso^J j. 



Fig-6(b): 



h 1 B± = f° dse ±icos e- is{EB i- EB)/n (37) 

J — OO 

= [° d se is (±^-E Bl +E B -iO+)/h 
J — oo 



i(±hco-E Bl +E B -iO+) 
h 



e is(±hw-E Bl +E B -io+)/h 



Fig.7(a): 



i(±fto; - E Bl + Eg — iO+) 



W = -E<b(]^) 1/2 ^] (38) 
/° /° dse ±iws 'e~ is ' {EA 2~ EB)/n+is{EA i~ EB)/n 

J — OO .7 — CO 

exp(- £{ J coth ^[(^) 2 + (^ 2 ) 2 + (^ - ^) 2 ] 

+ ^ ^ 2 - ^ 3 ) [coth cos s'o; a - i sin s'o^] 

+ ^a 3 (^ 2 - ^ 3 ) [coth ^ cos so;* + i sinsa;*] 

+ \0* 2 0t [coth cos(s - s + i sin(s - s')a;«] }) ■ 
Fig.7(b): 



I 



A 3 B 1 B± 



L<b( 



3 ° v M„/a; ff / y 2 



(39) 
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f° ds f° dse^^-E^/n-rs^-E^/n exp{ _l £ coth ^ [( ^ 2)2 + ( ^ 2 _^ 1)2 + ( ^ 1) 2 

J -co J —co 4 ^ Z 

exp(- 3 ) 2 [(N a + 1) + ^T** + ^±1***]) -I - [coth ^ coss'o,. - f si 



Fig.7(c) 



Ia 2 b 3 b± = f° ds f° dse ±icos ' 

J —CO J —CO 

E v od ( ft V 
B 3 sl M ,o; / j 2 



sins' cox] 

+ \ E^ 1 ^ 2 I coth ^ cos ( s + S 'K - »' sin(s + s')w K \ 

DC ^ 



(40) +^E^ 2 (^ 2 -^)Mh^ 

e is(E B3 -E B )/h-is'(E A2 -E B )/hy K oc' / ft a 1 /2 ^ 

ivl iys -0L IUJ 0L' 



— — cos stOcc — i sin sco a \ } 
(43) 



(z coth i — ^ — [sm(s — s jay — sin say J 
+ cos(s — s')ay — cos say) 



)/h 



— coth i — ^ — cos sa^ + z sin sa^J }, 



2 

The time integral of Fig. 7(d) is 0. 
Fig.8(a): 



Fig.8(c): 

Sib 2 a,b±(co,T) = f° dse-W**-** 

J —CO 

f° ds i e ±i(s+s')w e -i{s+s')(E Al -E B )/h 

J —CO 

exp{-iEcoth^(^) 2 



(44) 



Sia 2 a iB± KT) = /° dse ±- e --(E A2 -E^)/^ (41) Fig.8(d) 

J —00 

-00 



J— 00 



sins'a; a ]} 

J —CO 



(45) 



{(# -<, 2 )[/coth^sin S W -cossW] / ^-'^(^^^^[E^'b (— ^— ) 1/2 ^ 

J-00 ^ 2 1 M a /ay 2 

+<# [z coth ^ sin(s + s')ay - cos(s + s>v] - <# }] (cos ^ _ ■ CQth ^ sin ^ _ 1}] 
ex P {4 2 - ^) 2 + (^ 2 ) 2 + (^) 2 ] coth ^ 

exp{-im Al ) 2 coth^ 
+ 7 L^Qot (Qot - Qoc ) (coth — cos sa^ - zsmsa^J a 

— - 2^(^ ) (COth —— — COSS COpc — I Sins 0> a )}. 

+\ E&O* (coth ^ cos(s + s>« - i sin(s + s')cv a ) 

Fig.9(a): 

- 5 E « 1 " « 1 ) ( coth cos s '"« " ' sin s '"« ) > 
Fig.8(b): 

Sa 2 a 1 b±(^T)= /° dse"^"^ 

J —00 

J— o 



Sia 2 b iB± (^T)= /° dse ±^ e --(E^ 2 -E Bl )A 

J —CO 



)/h 



(42) 



/° d S ' e - i ( s+s, )( E »i- E B) / *{r< B (, 

J — CO • i 



? -i(s+s f )(E A -E B )/hr Ar / / T x fihcVct , f . 

e 1 Wa 2 AiB(S/S , i j, ^ zco th^ — sins co^i — coss co^) 

, -(zcoth^^sin(s + s / )o;^-cos(s + s / )a;^)]} 
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exp{-iEcoth^(^) 2 



/ ds'e ls c ° 1 cos s'oo, to = to 0L >, to\ = {E Bl — E B )/h 
J —00 



+ I E(^ 2 ) 2 (coth ^ cossa; a - i sinsa;,)} (46) /° dsV^ coss'a; -> /° dsV v 

2 ^ 2 .7 -co .7 — oo 



(^i-^'coss'a; 



Fig.9(b): 

S 2 A 2 B lB± (^T)= /° rf Se -<*(EA 2 -E Bl )/& (47) 

.7 —00 

/"° ds l e ±i(s+s')co e -i(s+s') (E Bl -E B )/h 

J — OO 



CO 



,A 2 



M a /ay 2 



, . , Bhco„i . 

(cos say — i coth ; — - — sin say — 1 )\ 



exp{-iE(^ 2 ) 2 coth^ 

z a z 
- ^X](^ 2 ) 2 [ cotn cos sc °a -i sin scope]}. 

The time integral of Fig.9(c) is 0. The time integral of 
Fig.9(d) is 0. 

3.3 Order K 2 contributions Fig.lO(a): 

rO rO rO , . „ 

Ib 3 b 2 b 1 b± = ds ds' ds"e ±l ^ (48) 

J —oo J —oo J —oo 

ex V {-is"{E Bl -E B )/h + is{E B3 -E B )/h + is'{E Bl -E B )/h} 



= / ds' cos s'(toi — z'0 + ) coss' 

J — oo 
rO 

+i / ds' sin s' [co\ — z'0 + ) cos s'to 

J — oo 

1 f° 

= - / ds' [cos s'(a;i — z'0 + +o;) +coss / (o;i — z'0 + — to)] 

2 J —CO 

+ - / ds^sins^o;! — z'0 + +co) +sins / (o;i — z'0 + — to)] 

2. J — oo 

1 1 

= — [sin s' (to i — z'0 + + to) 

2 to x - z0+ + to L v 1 7 

+z(-)coss / (o;i -z'0 + +o;)] 
1 



1 

+ - 



■[sins'^i — z'0 + — a;) 



2^1- z'0+ — to 

+z( — ) coss^a?! — z'0 + — a;)] 

: — }— [cos s' {co\ — z'0 + + to) + i sin s' (a?i — z'0 + + a;)] 

2 c^i — zlr + to 

1 1 

— [cos s'{co\ — z'0 + — to) +zsins / (a;i — z'0 + — to)] 

2 — z0 + — to 

_j gis'(coi—iO + +co) j gis'(coi—iO + —co) 

~ ~2 ' to 1 - z'0+ + to 2to x - z'0+ - to 

-l r 1 1 



+ ' 



(coth i — - — cos s a; a / + z sin s ay 

adiabatic introduce interaction 

E B3 - E B -> E B3 — E B — z'0 + 



consider 



2 to\ — z'0+ + to to\ — z'0+ — to 



'^-^sins'ov 



/•0 

/ dse !S 

J — OO 



(E B -E B )/h 



h 



i(E Bs -E B ) 



±hto — E Bl + Eg —> ±hto — E Bl +E B — z'0 + 
f° ds n e ±ius"-is"{E Bl -E B )/h = _ ft 

J-oo /( 

consider 



z'(±fto; - E B i + £ b) 



/ dsV V ( EB 2- £ B)/^ COSS / a ; a/ 
J —00 



J — 00 

rO . , rO 

/ ds^ 15 a;i sin s'to = / ds 7 cos s'toi sin s 7 a; 

J —oo J — oo 

+i I ds' sins^i sins'a; 

J — 00 

1 r 

= - / ds^sins^a;! — z'0 + +a;) — sins'^i — z'0 + — a;)] 

2 J —oo 

1 r° 

+ - / ds 7 [cos s r (6(7i — z'0 + — to) — coss'(toi — z'0 + +to)] 

2 J— 00 

coss^a?! — z'0 + + a;) i sin. s f (to\ — z'0 + + a;) 
2(to 1 -z'0+ +a;) 2(a;i - z'0+ +a;) 

^coss r (a;i — z'0 + — a;) z sins^a^i — z'0 + — to) 



2(to 1 - z*0+ - to) 2(to 1 - z*0+ - to) 
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gis'(a)\— iO + +cv) gis'(cvi—iO + —cv) 
+ 



2(co 1 - z'0+ + co) 2(cv 1 - z'0+ - to) 

lr 1 In 



2 <jO\ — z'0+ — a; a>i — z'0+ + a; 
in summary 

[° ds f e isf ^2- E B)/\ coih C ossW + zsins'ov 

J — 00 2 

— z /5/zav , 1 1 i 

~T COt 2 - z'0+ + a; + ~oo\ - z'0+ - cJ 

f r 1 1 
+ - 



2 — z'0+ — a; c^i — z'0+ + or 

Fig.lO(b): 

Ki B3 B 2 B 1 B±(^r)= /° ds e - /s ( EB 3- Efi 2) 

J —00 

/° ds'e- /(s+s ' )(EB 2- £ii i )/?i (49) 

J — 00 

Z* ds // e ±z , (s+s / +s // )a; e -z(s+s / +s // ) (£ Bl -E B ) A 

J —00 

'Z-r A B3B2 A B 2 Bi 



2V 



[coth i — ^ — cos s co a — i sin s co^. 



Fig.lO(c): 



dse- is{EB 3- EB z )/n 



f 

J — oy 

/*° ds / g ±i(s+s / )^ 6 - I '( s + s O ( £ b 2 -E Bl ) A 
J— oo 

f° ds " e -^s' + s")(E Bl -E B )/n (5()) 

J— 00 

2 L-s K B 3 B 2 — 



[coth ft^* cos(s' + s")a; a — z sin(s / + s")a; a ]. 
Fig.lO(d): 

J —00 

/° ds ' e -'(s + s')(E B2 -E Bl )/h (51) 

J —00 

/° cfs"e" ! ' (s+s ' +s " )(EB i" EB)A 

J —00 

9 Z-f A B 2 Bi A BiB" 



2 ^ w 

r i-u /^^a // • • // i 

[coth — - — cos s a; a — z sin s a; a J . 



4 Asymptotic calculation of the time integrals 

To illustrate how to carry out time integrals, we take 
Eq. ^19)) as an example. Extend the integrand to the 
whole complex— s plane and write s = u + iv. In 
(u,v) language, 

coth cos o;s + z sin o;s 

= cos (coth coshow — sinhow) 

+ i sin oju (cosh a;i7 — coth sinh a;i?) . (52) 

At one end point (0, 0) of the integral, the imaginary 
part of Eq. (|52)) is zero. There are two paths on which 
the imaginary part of Eq. (|52)) is also zero: (1) u = 
and (2) v = fth/2. They are the steepest descent paths. 

Because the integrand is analytic in the whole 
complex— s plane, according to the Cauchy theorem, 
we can deform the original contour C: (— oo, ())—>• 
(0,0) to steepest descent path Q+C2+C3, where Q: 
(-X,0) (-X,j8ft/2),C 2 : (-X,ph/2) (0,j8ft/2), 
C 3 : (0,]8ft/2) (0,0), In the end, let X 00. Obvi- 
ously, } c = J Ci + / C2 + / C3 . 

Along path C\: s = — X + iv, where v: —> fth/2. 
The integral becomes 



-iX(co—co f A 



Jo 2^ 



[cos a> a X(coth ^^ c ' g cosh a> a z> — sinha; a z>) 

+ ism avX(coth sinha> a p — cosha> a p)]}. 

(53) 

Because X is a large positive number, the phase of the 
integrand changes sign wildly due to sino^X, and 
Eq.(|53| is negligible. 

Along C2, s = u + z^, where u: —X — » 0. Ex- 
cepting the factor exp{- £ a ( ^~^ )2 coth the 
integral becomes 



, phw/2 e Pfa>' AiA /2 



f 

J 1 



due t«u-iuu AiA 



ex P {^E(^ -^) 2 cscfc^cosa, aM }. (54) 
The main contribution comes from the neighborhood 

2 

of u = 0. Expand cos o^tz « 1 — w 2 ^, then Eq. (|54)l is 
proportional to 



j icvu—iucv' A A —Cu 2 

due A i A e LM , 

-00 



(55) 
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where 

Around u = 0, the Taylor series of e lwu lUL ° A \ A con- 
verges rapidly. One can then easily carry out (|55)) . 

Along C^: (0,0) —> (0, the integral becomes 

i T 2 dve-*"-"'^ exp{i - tf) 2 

Jo 2^ 

(coth ffi^ cosha^z? — sinho;^)}- (56) 

According to the Laplace method, the main contribu- 
tion comes from the neighborhood of v = ph/2. The 
leading behavior of (|56)l is then 

exp{£^-^) 2 csc^}. (57) 



5 Time integrals for Hall mobility The time inte- 
gral for Fig. 11 (a) is 



,C 

Uaa 2 a 1 — / 

.7 — i 



dse' 



is(E A -E A 



(58) 



1 P ht0 « 



[ dsV V( ^- £ ^ ) exp(^{-jCoth 

J-oo a 4 

[(^ 2 - tf 1 ) 2 + (^ 2 - ^) 2 + (^ - e£f] 



+^{et-et){et-et) 

[COth — - — COS SCO a — 1 Sin S(V a \ 

[cos SCO a — coth - — - — i sin sco a \ 



[cos(s + s')o; a — coth ^^zsin(s + s')a; a ]}). 
The time integral for Fig. 11(b) is 

U' AAlA2 = f° dse is ^- E ^ /n (59) 

J —00 



^-^-Ea)/* exp(0" | coth 

[(^ 2 - + (et 2 - e^) 2 + (^ Al - e^f] 
+\{^ 2 -^){et 2 -et 1 ) 

[coth i — ^ — cos so^ + i sin sa; a J 
[coth ft^* cos(s + s')^ — z sin(s + s')a; a ] 

[coth i — ^ — coss a; a - /sins a^J jj. 

The topology of Fig. 12(a), Fig. 13(a) and Fig. 12(b) 
are the same. They have same time integral: 

Z AA3A2Al = f° dse- is ^- E ^ /n (60) 

J — CO 

f° ds'e- is '^i- E ^ /h 

J — oo 

f ds " e -"^ 3 -E A )/n exp( £ { _ 1 coth P?*L 

[(0a 1 - ^ A ) 2 + (e« 2 - o^) 2 + (^ - et) 2 + (^ 3 - 0*) 1 ] 

-l(9a 2 -0^)(e^-e^) [coth ^ cos o; a s' - i sin a> K s'] 
-\(e« 3 -0a)(O« 2 - a" ) [coth ^ cos a> a s" + i sin u> a s"} 



+ ^(0a 2 -0a 3 )(0a 1 -0^) 

[coth coso; a (s + s f ) — ismcv a (s + s f )] 
+ ^ 2 - 0a 3 ) (0« 2 - 0« 1 ) [coth PUp. cos cv K s - i sin o; a s] 



-0$) 



[coth J — ^ — coso; a (s + s — s ) — i sina^s + s — s )\ 



+ \{0a 3 -0*)(0a 2 -0a r ) 

[coth ^ l ^ DL cosa; a (s — s") — ismco^s — s")]}). 
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The topology of Fig. 12(c) and Fig. 12(d) are the 
same. They have same time integral: 

Z' AA3A2Al = f dse is ^- E ^ h (61) 

J —00 

f° ds'e is '^- E ^ /n 

J — CO 

f° ds"e- is "^- E ^ /h 

J — 00 

exp(0-Jcoth^[(^-^) 2 

DC 

1 A 2 _ a A l\2 , f P A 2 _ nA 3 ^2 , / P A 3 _ n A\2] 



+m 2 - cr + (^ 2 - + (c 3 - m 

(a A 2 q a i\(qA\ qA\ 



[coth — - — cos CjOqcS + z sin co^s \ 



1 



r -J ShCOtX f, . . //-, 

[coth i — ^ — cos o^s — z sm o; a s J 

+\{& tit) 

[coth ffi^ff coso^(s + s') + z'sina^s + s')] 

\ (Q A 2 a A 3\(r\A 2 

[coth i — ^ — cos a^s + z sm a; a sj 

+\{&-ti){f£ 1 -tf) 

[coth fi^ 7 * cosa^(s + s' — s") + isina; a (s + s' — s")] 



[coth ^ ~ cosa; a (s — s") + z sina; a (s — s")]}). 

The topology of Fig.l3(b), Fig.l3(c) and Fig.l3(d) 
are the same. They have same time integral: 

Y A3 a 2 a iA = f dse- is ^- E ^ ,n (62) 

J — 00 

[° d s > e -i(s+s')(E A2 -E Al )/h 

J — CO 

f ds " e -^ + s")(E Al -E A )/H exp{ _ 1 EcQth 
J— oo 4 ^ 2 



[coth ffi^ cos(s / + s")cOcc — i sin(s / + s")cOa\ 

+ J E(^ 2 - ^ 3 ) (^ 2 - 0* 1 ) I coth ^ cos - * sin 



i 



2~ 



[coth cos(s + s f + s")^ — sin(s + s f + s")^] 



[coth ^^_r cos(s + s')a; a — z sin(s + s')co,x\ 
- \ E(0« 2 - & ) 3 - # ) [co* cos scv a - i sin sco a ] }. 



6 Potential energy partition For extended states, 
we separated h e into two parts: h e = h® + h e _ p j lr 
where 



^2 A/" 

■E 

n=l 



2^V?+E^(r-^n), (63) 



and 



A/" 



3Af 



dU 



^ = -[u n -V r U(r-K n ) = (64) 

n=l y=l 0A y 

The eigenvectors £g and eigenvalues £g of h® are 
taken as those of the Kohn-Sham Hamiltonian. To 
make partition (|63I64|I work, one needs: (1) h e _ p j l does 
not change the gross feature of the phonon spectrum 
and the normal modes; (2) the changes in the eigen- 
values and eigenvectors of h® caused by h e _ p j l are 
small. 

For a localized state, we partition the potential en- 
ergy as 

E^-^<)= E U{r-n n ,u*) 

n nGD Al 



+ £ U(r-1l v ,vP v ). 



(65) 
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Because the wave function of a localized state is con- 
fined in a finite region, the e-ph interaction from the 
atoms outside D A is negligible: 

£ U{r-K v ,u v v )^ £ U(r-n v ). (66) 

In partition (|65)) , the zero-order Hamiltonian for lo- 
calized state A\ is 



h» = 



1 2m 



V?+ £ U(r-Kn,u*). (67) 

nGD Al 



Denote the lowest eigenvalue and corresponding 
eigenvector of h° A as E Al and (p Al , using stationary 
perturbation theory, we have (p Al ~ <p\ , and 

^ = 4," E dpASrl,, (68) 

where and E® A are the Kohn-Sham eigenvector 
and eigenvalue. 

d PAi = Jdv\<p° Ai (r)\ 2 dU/dX pAi , (69) 

is the e-ph coupling for the p A vibrational degree 
of freedom. Here dU/dXp A ^ is the derivative of the 
single-electron potential energy U respect to the 
nuclear degree of freedom. 

The partition ([65]) is justified if 

| ( drf Ai {r-n Al ) £ U(r-K p )(p Al (r-K Al )\ 

<l E d VA^VA x l (70) 
Pa^D Ai 

and 

| fdr^t) E U(r-K p )cp Al (r-K Al )\ 

<l E d PA^p Al l ( 71 ) 
PA r £D Al 

One can easily estimate the RHS of Eqs. (|70l71[) , the 
total e-ph interaction energy for a localized state (p Al . 

The number of atoms inside D Al is (£/d) 3 , where £ 
is the localization length of ^ . The RHS of ([70)) is 

= (47re ) _1 (Z*e 2 /d)(£/d)(^/d). The attractive en- 
ergy between a nucleus in D Al and the carrier in D Al 



is ~ (4zreo) 1 Z*e 2 -/R\2, where R12 is the distance be- 
tween the centers of two localized states. The LHS of 
dZQ|is 

^ 2 /dfe- R ^ /d (4:7te )- 1 Z^e 2 /R 1 2 / so that ££8 is sat- 
isfied. The LHS of EQ) is 

(£/L) 3/2 (£/d)(4zre )- 1 Z*e 2 /d, where Z* is typical 
effective nuclear charge, L is the linear size of sample. 
For a well-localized state £ << L, and factor (£/L) 3 / 2 
assures that condition (Q} is satisfied. Since parti- 
tion ([65)) includes polaron as an example, we can also 
apply the present ansatz to obtain the polaron states 
from first principles by computing the static displace- 
ments caused by e-ph interaction. 

7 Scenario for simulation We describe how to 
implement the present method in ab inito code. For 
an amorphous semiconductor, using a super-cell 
structural model, one first calculates the eigenvectors 
{<Pq} and eigenvalues {E^} of Kohn-Sham Hamilto- 
nian h e at r point. From the vibrational Hamiltonian 
h v or dynamical matrix, one can obtain vibrational 
spectrum {00^} and eigenvectors {A moc }. 

For a given state (p® c , the electron-vibration cou- 
pling constant for the p th nuclear degree of freedom 
is given by: 



d c v 



j dr\(p c (r)\ 2 dU/dX p , p = 1,2, • • -3AT, (72) 



where U is the effective single-particle potential en- 
ergy. The static displacement x m of the tn th degree of 
freedom is 



p 



)mpj 



1,2, •• -3Af, (73) 



where k 1 is the inverse matrix of the dynamical ma- 
trix k. The vibrational amplitude x m of the m th degree 
of freedom is determined by the amplitudes {©#} of 
normal modes: 



x m = E A ^ *' ®<* max{ 1 



' kT 
Mo? 2 ' 



h 



MiOoc 



h 

(74) 

where A moi is the minor of the determinant 
\Kk — ^MiSj^ = 0. If for some indexes x m > x m , we 
say state (p^ is a localized state (p A . If none of the three 
degrees of freedom of an atom satisfies x m > x m , we 
say that this atom does not belong to D A . Otherwise 
we say that the atom belongs to D A . We will treat 
the e-ph interaction in zeroth order together with H v 
and consider the polarization effect expressed by x^ m . 
The polarization effect of a localized state is like a 
small polaron, but the electronic localization in AS is 
caused by the static disorder in h® rather than h e _ p j l . 



Copyright line will be provided by the publisher 



pss header will be provided by the publisher 



15 



The topological disorder contained in {1Z n } is much 
larger than the static displacement {x%} induced 
by e-ph interaction, the binding force caused by the 
static disorder is much larger than the self-trapping 
force caused by the e-ph inter actiontl]. If for a state 
(pc, none of vibrational degrees of freedom satisfies 
*m > % m, we say this state (p^ is an extended state 
and treat e-ph interaction as perturbation. 

From {x® A }, the shift a 1 of origin for the oc th 
mode is then obtained [1 J. With the shift in origin 

{O^ 1 } for each localized state, the time integrals in 
Sec|3| and [5] can be found numerically The small 
parameters /, J 7 , K' and K for the residual interac- 
tions are computed Q by the expressions below [text 
Eq.(16)]. Finally one obtains conductivity in tables 
4 and 5 and Hall mobility in [text Eq.(23) and table 
8]. The present results are applicable to crystalline 
materials and amorphous metals if the scattering 
mechanisms are only the disorder potential and the 
electron-phonon interaction. Although we did not 
explicitly deal with the scattering effect of the static 
disorder, the scattering of Bloch waves by the disor- 
der potential is already taken into account by [text, 
Eq.(17)]. The reason is: the basis set for xp(°\ 

xp^ and xp^ are solved from the Kohn-Sham Hamil- 
tonian, so that xp^\xp^ and xp^ are fully dressed by 
the disorder potential. 
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In this paper, we review and substantially develop 
the recently proposed "Microscopic Response Method", 
which has been devised to compute transport coeffi- 
cients and especially associated temperature dependence 
in complex materials. The conductivity and Hall mobil- 
ity of amorphous semiconductors and semiconducting 
polymers are systematically derived, and shown to be 
more practical than the Kubo formalism. The effect of 
a quantized lattice (phonons) on transport coefficients is 
fully included and then integrated out, providing the pri- 
mary 



temperature dependence for the transport coefficients. 
For higher-order processes, using a diagrammatic expan- 
sion, one can consistently include all important contribu- 
tions to a given order and directly write out the expres- 
sions of transport coefficients for various processes. 
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1 Introduction The career of O. F. Sankey exempli- 
fies the ideals of scientific exploration, and even adventure. 
His many contributions to semiconductor, materials and 
biophysics are testimony to a man of ingenuity, energy and 
integrity. He has also been a rigorous but patient mentor to 
many contributors to this volume, including D.A.D. This 
paper on transport in complex materials is offered with af- 
fection and gratitude on the occasion of his sixtieth birth- 
day. 

The Kubo formula^ has been used to calculate the con- 
ductivity and Hall mobility of small polarons in molecu- 
lar crystals ll5ll6ll7ll8l [MT0ll . The results obtained by differ- 
ent authors are inconsistent. The reasons are: (i) the imagi- 



1 In this paper the Kubo formula refers to the original treat- 
ment of Kubo[T), not the subsequent approximate form of 
Green woodlEl. which is discussed elsewhere(3p1l. 



nary time integral in the Kubo formula is complicated; and 
(ii) there is no systematic way to classify various transport 
processes induced by the external field and by the resid- 
ual interactions. The small polaron model is inadequate 
for amorphous semiconductors (AS) in two aspects: (i) the 
states at the fringes of valence and conduction bands are lo- 
calized by topological disorder. At moderate temperature 
both localized and extended states are accessible to elec- 
tronic transport processes iTTA . Fig[T]is a schematic energy 
spectrum of AS; (ii) A carrier in a localized state (j) A po- 
larizes the nearby atoms: the static displacements d PA of 
the atoms (pa is the index of atoms participating in local- 
ized state (j) A ) induced by the e-ph interaction is compara- 
ble to or even larger than a v ,the amplitude of thermal and 
zero point vibrations. In contrast, the static displacements 
of atoms caused by the carriers in extended states ^ B are 
negligible O. 
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Figure 1 Schematic energy spectrum of amorphous semi- 
conductors or semiconducting polymers. 



The imaginary time integral in the Kubo formula can be 
avoided in the microscopic response method (MRM)|[T3ll3l 
4 ] . If an external perturbation on a system can be expressed 
via additional terms in the Hamiltonian (a "mechanical 
perturbation" m, electromagnetic field is an example), the 
microscopic response can be obtained from the continuity 
equation for charge, the time-dependent Schrodinger equa- 
tion and the initial conditions. The measured macroscopic 
response is a coarse-grained average, and an ensemble av- 
erageof the microscopic response! 13]. Thus one does not 
need to calculate the macroscopic response by averaging 
the response operator over density matrix, and the imagi- 
nary time integral is avoidedO. The purpose of this paper 
is to review the MRM and to report the full formulae for 
the conductivity and Hall mobility. 

The MRM is equivalent to and simpler than the Kubo 
formula[3|. Because the microscopic current density in a 
state is a bilinear form of the wave function of that state 
connected by the current operator, each transport process 
is a temporal evolution of a state under the driving force 
of the external field and the residual interactions which are 
time ordered in a specific way. The contribution to the mi- 
croscopic current density from a process is a product of 
several transition amplitudes and a connector of the cor- 
responding current operator. One may visualize the trans- 
port processes by a series of diagrams. The topology of 
diagrams gives us a systematic classification about various 
transport processes. It becomes easier to compute transport 
coefficients to a given order of residual interactions. 

The consequence of strong polarization was first dis- 
cussed by Marcus for electron transfer in polar solvents lfT4l 
and by Holstein for small polarons in molecular crystals [15) . 
Emin and others |[T6l have used small polaron theory to 
compute transport properties for well-localized carriers. 
However two related features in this model do not agree 
with the nature of amorphous semiconductors: (i) elec- 
tronic localization in an amorphous semiconductor is 



caused by the static positional disorder rather than the 
strong e-ph interaction; (ii) the localized states caused by 
the geometric disorder involve several atoms rather than 
one atom lfTTl . 

Miller and Abrahams studied the carrier hopping in a 
doped crystalline semiconductor. They assumed that the 
carriers are trapped in impurity or defect states close to 
the Fermi level, and the e-ph interaction is weak.Then the 
transitions between donor sites is mainly affected by the 
single-phonon absorption or emission lfT8l . 

Traditionally transport theory takes the zero order 
Hamiltonian Ho as a sum of the vibrational part H v and the 
electronic part H e in the force field of fixed nuclei, and the 
e-ph interaction H e - V h is viewed as a small perturbation: 



^0 — He + Hvib, H — Hq + H e 



-ph- 



(1) 



We will use h e and h e - p h to denote the corresponding 
single-particle quantities. If the e-ph interaction is so weak 
that the static displacements of atoms induced by the e- 
ph interaction is negligible compared to a v , partition (QJ is 
reasonable. Both localized states and extended states are 
eigenstates of h e , H e - V h is the unique residual interac- 
tion and induces [19 1 all four possible transitions among the 
eigenstates of h e : (i) from a localized state to another lo- 
calized state (LL), (ii) from a localized state to an extended 
state (LE), (iii) from an extended state to a localized state 
(EL) and (iv) from an extended state to another extended 
state (EE). 

Both experiments and ah initio simulations show that 
in AS the e-ph interaction for the well-localized tail states 
is strong (d PA > a v ), while the e-ph interaction for the less 
localized states and for the extended states is weak ||20l 
ED • One has to adopt different Hamiltonian partitions for 
the two situations. For a well-localized state A (we use A 
with or without subscript to label a localized state), we take 
Hq g = K + Va as the single-particle Hamiltonian, where 
Va is the effective potential energy of the ions inside the re- 
gion Da where (j) A (eigenfunction of h e ) is nonzero. One 
can easily show that (j) A « (jP A even under perturbation 
h e - p h, where (j) A is the ground state of h$ e lfl2ll . The zero 
order Hamiltonian is ho = h$ e + (H v + h e - p h). The at- 
traction from other nuclei outside D A is taken as a pertur- 
bation. In this partition, LL and LE transitions are caused 
by the transfer integrals 1 12], cf. Sec. 6 of Ref. ll22l . For an 
extended state £ B (we use B with or without subscript to 
label an extended state), the zero order Hamiltonian was 
taken as hs = h e + H v , h e - p h is the perturbation. EL and 
EE transitions are caused by the e-ph interactions lTl2l . 

In this paper, we apply the MRM to compute the con- 
ductivity and Hall mobility in amorphous semiconductors. 
In Sec 12 the measured macroscopic response is obtained 
by taking a spatial and ensemble average over the micro- 
scopic current density. The required input is the many-body 
wave function \P' (t) of N electron + J\f nuclei in an external 
field. By means of the single-electron approximation and 
the harmonic approximation for vibration, in Sec 13 one 
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can use a simplified many-body wave function ip'(t) to ex- 
press the spatial averaged current density, ijj'(t) describes 
the motion of single electron in coupling with the nuclear 
vibrations. Expanding ip'(t) with localized and extended 
states, the evolution equations of the transition amplitudes 
can be derived from the time-dependent Schrodinger equa- 
tion. It is convenient to compute the transition amplitudes 
in normal coordinate representation. In Sec|4l we outline 
how to obtain ip\t) to required order in residual interac- 
tions and in external field. In Sec|5j we compute the con- 
ductivity from both localized and extended states. With 
the help of a systematic diagrammatic perturbation expan- 
sion, one can determine conductivity to any order in ex- 
ternal field and small parameters of residual interactions. 
We point out why some important contributions have been 
missed in previous calculations of small polarons based on 
the Kubo formula. The non-diagonal conductivity needed 
for Hall mobility is calculated in Sec|6l To describe the 
Hall effect, one has to expand the current density to the 
2 nd order in the external field: one order is electric field, 
another is magnetic field. The MRM shows that: (i) there 
is a new type of term in the current density; (2) there is 
an intrinsic interference effect between electric field and 
magnetic field. A new type temperature dependence is pre- 
dicted. In the MRM, the atomic vibrations are described 
by quantum mechanics, the results obtained are correct for 
any temperature. 

2 Macroscopic current density Consider an amor- 
phous semiconductor, with N electrons and M nuclei in an 
electromagnetic field described by potentials (A, 0). De- 
note the coordinates of the N electrons as ri , T2, ■ • • , r/v, 
the coordinates of J\f nuclei as Wi, W2, • • • , Wjv". The 
state \P'(t) of the system is determined by the Schrodinger 
equation 

ififfl'/dt = H'V, (2) 

where H' = H + Hf m is the total Hamiltonian, H is the 
Hamiltonian of system, Hf m is the field-matter interaction. 
The arguments of \P' are (ri, • • • , r/v; Wi, • • • , W;v; t). 
In this work, we focus on the transport coefficients far be- 
low the frequency of infrared radiation. The direct contri- 
butions from nuclei will not be written out. For an ac field 
with higher frequency, especially in the infrared range, one 
has to take into account the direct contribution of the mo- 
tion of atoms. 

The microscopic current density in state \P'(t) 
isOH 

= --A(r;t)n'(r;t) (3) 
rn 

2m J 

where n f (r;t) = N f dr'^'*^' is the carrier density in 
state \P'(t), dr' = dr 2 • • • dr N dWi • • • dW^, the argu- 
ments of \P' in Eq.® are (r, r2, • • • , r/v; Wi • • • Wjv"; t). 
The first term of Eq.© expresses the oscillation of 'free 



electrons'. Comparing to the second term, its contribution 
to transport coefficients is negligible (23]]. We will not dis- 
cuss it further. 

By averaging over a 'physical infinitesimal' volume 
element (HO i? s around point s, the coarse-grained cur- 
rent density is : 

if (s,t) = Q- 1 [ drjf (r,t), a = x,y,z. (4) 

For a mechanical perturbation, state \P'(t) is determined by 
the initial condition H/'(—oo). Because one does not know 
what state the system was initially in, one must average 
j% over all possible oo) to compute the measured 
macroscopic current density j a : 

j a (s,t) = J2 p [*'(-oo)}jt (s,t), (5) 

where P\&'{—oo)\ is the probability that the system is 
in state \P'(—oo) before the external field is adiabatically 
introduced. P[\P f (— oo)] depends only on the energy of 
state \P'(—oo) 9 and may be taken to be the canonical 
distribution!!!!!. 

3 Harmonic approximation and single-electron 
approximation In the solid state, W n = 7l n + u n , 
where^ n = (X 3(n _ 1)+1 , X 3(n _ 1)+2 , X 3n ) and 
u n = (x 3 ( n _i) + i,X3( n _ 1 )+2,^3n) are the equilibrium po- 
sition vector and the vibrational displacement vector of the 
n th nuc i eus respectively. In the harmonic approximation, 
the vibrational Hamiltonian is 

H * = E -^M,h + \ £ kjkXjXk > (6) 

3 3 jk 

where (kij) is the matrix of force constants. 

Because in AS the correlation between electrons is 
weak, one can use the single-electron approximation to 
(t) . The arguments of the simplified single electron wave 
function ip f (t) include only the single electronic coordinate 
and the vibrational coordinates of nuclei. The state ?// of a 
carrier in an external field satisfies: 

iHdil> , /dt = tiip , (r,{x j },t), (7) 

where r is the coordinate of the carrier, h! = h + hf m 
is the Hamiltonian of [system + external field]. In gauge 
V- A(r,t) = 0, 

hf m = (^e/m)A(r)-V r + e 2 A 2 (r)/(2m) + e0(r), (8) 

is the coupling between the carrier and the external field. 
Here h = h e + H v is the Hamiltonian without field, where 

-fi 2 M 
h e = —V 2 + Y / U(r-Tl n - u n ), (9) 

n=l 
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is the electronic Hamiltonian. U(r — 7Z n — u n ) is the effec- 
tive potential energy between the electron at r and the n th 
nucleus. Denote the 3AT vibrational degrees of freedom as 
{xj, j = 1, 2, • • • , 3Af}, then the spatially averaged cur- 
rent density © at s is reduced to II 13 11241 



iheN e 
2mi? s 



/ dr [\\dxj] 



(10) 



4 Evolution equation in external field The prime 
ingredient required for the conductivity and Hall mobility 
is the single particle state ij)'(t) of the system in an external 
field |[T3l . To find ip' (t) , we expand it using the approximate 
eigenstates of the single-particle Hamiltonian lfT2l : 

^'{t) = Y j aA{t)4> A + Y J b B {t)t B , (11) 



where the arguments of ip' are (r; x±, X2 : • • • , xsj\f), a>A{t) 
is the probability amplitude at time t that the carrier is in 
a localized state (p A while the vibrational displacements of 
J\f nuclei are X2, • • • , x^. Here, b B (t) is the proba- 
bility amplitude at time t that the carrier is in an extended 
state £ B while the vibrational displacements of M nuclei 
are x\ , x^, ■ • ■ , x$m- If me external field is not too strong, 
it is convenient to put the change in ip' caused by external 
field into the probability amplitudes rather than in the zero- 
order eigenstates. Substituting Eq.(Qj} into Eq.©, the time 
variation of probability amplitudes under the external field 
are determined by fT2l : 



d_ 

dt 



and 



d 



[ih- - h Bl ]b Bl = J'b>a + £ K B 0t lB bB, (13) 

A B 

where 

jtot _ j _j_ j field ry-ftot _ rsl , ts-i field siA\ 
J A 1 A — J A 1 A^J AlA , J^A 1 B — JX A 1 B^ JX A 1 B ' v 1 ^ 

and 

J'^A = J' Bl A + 4:i\ K% B = K BlB +Ki% d , (15) 
and 

Ha 1 = H v + E Al , h Bl = H v + E Bl , (16) 

where transfer integral J Ai a = J drcp* Ai J2 p ^d a U ( r ~ 
7Z p )(j) A causes a transition from localized state <p A 
to localized state (j) Ai , the transfer integral J' BlA = 
J dr£* Bl ^2 p ^ Da U(r — Tl-p)4>A induces a transition from 



localized state (j) A to an extended state £ Bl . The e-ph in- 
teraction K AiB = J drcj) Al J2j x j(9U/dxj)^ B causes a 
transition from extended state ^ B to localized state <j> A . 
The e-ph interaction K BlB = J dr£* Bl J2j x j(9U/dxj)^ B 
induces a transition from extended state ^ B to another ex- 
tended state £ Bl . J A f A = j dr(j)* Ai h fm<t>A I s me coupling 
between two localized states (j) A and <j> A caused by exter- 
nal field, etc. 

The microscopic response is expressed by the state 
ip'(t) of [system + external field]. We replace the second 
line of Eq.flS with (^V r ^* - f*V r f) ~ (^V r V>* - 
^*V r ?/ ; ), where ^(r, {xj} : t) satisfies ihdip/dt = hip, 
is the state of system without external field. (ipV Y ip* — 
ip*V r ip) represents the microscopic current density when 
no external field is applied to the system. For carriers in lo- 
calized states and in extended states (scattered by static dis- 
order), the contribution to the coarse-grained current den- 
sity from (ipV r ip* — ip* Vr^) is zero. The spatially av- 
eraged microscopic current density to second order in the 
field is[ 13] 



N e he 



/ dr [[[dxj 



(17) 



{Im(V (0) V r V (1) * - ?/> (1) *V r ?/> (0) ) 
+ Im(V> (0) W> (2) * - ^ (2) *W (0) + V (1) V^ (1) *)}, 

where ip^ is the change in state to first order in field, and 
ip^ is the change in state to second order in field, ip^ 
is the state of carrier at time t without an external field. 
Hereafter, the superscripts (0), (1) and (2) on ip indicate 
the order of external field, the order of residual interactions 
will be denoted as subscripts. For example ip ^ represents 
the change in state to first order of external field and to first 
order of J. If not explicitly stated, we understand that ip^ , 
ip^ and ip^ are fully dressed by the residual interactions 
J, J', K' and K. 

Unlike the nearly free carriers in crystalline materi- 
als or the extended states of an amorphous semiconduc- 
tor, the carriers in localized states cannot be accelerated 
by an external electric field. Because the force produced 
by the external field is much weaker than the binding force 
from static disorder, the quantum tunneling probability l25l 

r = exp{-4(2m) 1 / 2 4 /2 /(3^e£;)} produced by an ex- 
ternal field is negligible, where £5 is the binding energy 
produced by static disorder. In addition, for a field of 10 5 
V-cm -1 and a typical distance 10A between two neighbor- 
ing localized states, the upper limit of the transition mo- 
ment is 0.0 leV, the same order as J and J' . Denote S as 
one of the small parameters J, J 7 , K' ', K from the resid- 
ual interactions. One must calculate ip^°\t), ip^\t) and 
ip^ (t) to order S° and S 1 and keep all possible order 5° 
and S l contributions in current density j. 
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To find ip'(t), it is convenient to transform the vi- 
brational displacements {x} of the atoms into normal 
coordinates {0} lfT2l . We expand probability amplitude 
a Al (■ • • O a • • • ; t) with the eigenf unctions of h Al ' 

£ C^(...N' a ... ] t)^e- i < L} /\ 

• • • • • • 

(18) 

where C Al (• • • N' a • • • ; t) is the probability amplitude at 
moment t that the electron is in localized state (j) Al while 
the vibrational state of the nuclei is characterized by oc- 
cupation number {N^a = 1, 2, • • • , 3AT} in each mode. 
Similarly, we expand &b x (• • • a • • • ; t) with the eigen- 
functions of Iib 1 : 

...AT;... 

(19) 

Eqs. (ll2ll3l) become the evolution equations for C A and 
F B . Using perturbation theory, one can find the probability 
amplitudes C A and F B to any order of field for any initial 
conditions. Substitute C A and F B into Eqs. (ll8ll9lh and 
one obtains the probability amplitudes a A and 6# in the 
normal coordinate representation. Substituting a a and bs 
back to Eq.flj}, one determines the state ?// (t) of system in 
an external field. The results are listed in Sec.l of Ref. l22l . 

Table 1 The origin of various conduction processes 



(^ (0) V r ^ (1) *-^ (1) *V r ^ (0) ) links all components of ^ (0) 
to all components of i^^* . 




a) b) 

Figure 2 Lowest order processes contributing to conduc- 
tivity. Initial state is localized state (j) A \ order J° and [J f ]° 
contributions to conductivity (see text). 




w w w w 

a) b) c) d) 

Figure 3 Initial state is localized state <p A \ order J 1 and 
order [J'} 1 contributions to conductivity. 



order 


expression 




diagrams 


s° 




-v4V*v r ^) 


Figiaia 


s 






FigEJEl 


s 






Fi2l4l5l8l9l 


K 2 




-^V*v r ^) 


FiglTOal 


K 2 






FisllObllOcllOdl 



5 Conductivity For the conductivity, one needs only 
the first term in Eq.JTTl). To order S° and S 1 , all possible 
contributions to j are classified in Table[T] If the initial state 
is a localized state <fr A , the expression for j can be found by 
substituting the corresponding if)^ (t) and ^ x \t) into the 
first term in Eq.JTTl). There are 14 processes contributing to 
conductivity. This can be understood as following. There 
are three terms in (1) the free evolution of state (j) A 
without external field and residual interactions; (2) J takes 
(\) A to another localized state; and (3) J 1 takes localized 
state (j) A to an extended state. As indicated in Eqs. (ll2ll3l) , 
e-ph interactions K and K' do not couple a localized state 
to any other state. There are two order S° terms in ip^\ 
the external field can bring <p A to either another localized 
state or an extended state. There are 8 order S 1 terms in 
they differ in the time ordering of the residual in- 
teraction S and external field. The current density operator 




a) b) c) d) 



Figure 4 Initial state is localized state <j) A \ order [J'} 1 and 
order J 1 contributions to conductivity 




a) b) c) d) 



Figure 5 Initial state is localized state (j) A \ order [J 7 ] 1 , K r 
and K contributions to conductivity 

We may use diagrams to visualize the 14 processes 
(cf. FiglJIQ. The expression for the current density 
Eqs. (11715b and the perturbative solution of Eqs. (ll2ll3l) 
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Table 2 States coupled by the current operator 



symbol 


expression 


A a/sa A ^ 




A->Bi 


/„. *(* A v rf - ttl v,^){4: ; Vi~° V'< i, - £r " ,)/ " 


B-^ A 2 


/„.({ B v r « 2 - «,v,f B )(a>« , l4"- , ) e "< ;:, - 4 "° ) ' / " 


B / w Bi 


«.-' /„. *« B v r G, - C V4 B )<4f ;, |s<"->>e" < < ; - , - 4 ''° ,)/ " 


A 2 — > A 3 





suggest some rules to write out the contribution to j for 
any process. ip^°\t) and ip^\t) are two ingredients of 
j. In Eqs. dl 11181191) , ip {0) (t) and ip (1) (t) are expanded 
with the eigenfunctions of h$ e + Ha and of ftoe + h<B- 
The expansion coefficients are transition amplitudes 
from a given initial state. The current density operator 
(^ (0) V r ^ (1) * - ^ (1) *V r ^ (0) ) links one component of 
ip^ to one component of ip^* (connector). Thus j is a 
sum of many terms, each term is a product of a connec- 
tor and one or several transition amplitude(s). We draw a 
wavy line from one component of ip^ to one component 
of xp^*. The spatial average produces Q~ x f reQ dr be- 
fore the microscopic current density, cf. the second to the 
fifth lines of Table [2 Because ip^ appears in j, the free 
evolution ip^}(t) of the initial state or one of two final 

components in ip^}(t) will appear in j. We use a solid 
arrow line upward (top right or top left) to express the 
transition amplitude in ^} (t) caused by one of residual 
interactions: the arrow points from initial state to final 

CD* 

state. For a process in which ip K s o appears in j, the com- 
plex conjugate (cc) of the transition amplitude(s) in ip^} (t) 
will appear. We draw a solid arrow line downward (lower 
right or lower left) to express the cc caused by an external 
field: the arrow points from the initial state to the final 
state, see Table [3] For a process in which tp^i* appears in 

j, the cc of two transition amplitudes in ip^} (t) will appear. 
We draw two successive solid arrow lines downward, one 
represents the cc of the transition amplitudes caused by 
an external field, another expresses the cc of the transition 
amplitudes caused by a residual interaction, see Table [3] 
By retaining only order S° and S 1 terms in j, FigEJ- \5\ 
describes all possible combinations. The two processes in 
Fig Gil result from ip^?o and tp^o *. The four processes in 
FigOresult from xp^} and ip^o* . The 8 processes in Fig.[4| 
and \5\ result from ip^o and *• 

With the help of the diagrammatic rules listed in tables 
[2|and[3l one can easily write out the corresponding macro- 
scopic current density for each conduction processes. Tak- 



ing Fig|3a|as example, the contribution to j is: 

(-^r) E /wnw E [wo^j] 

A—N a — oi A 1 ---N' (X --- 

£ [l-f{E As )\Q- 1 f dr(^ 1 V r ^,-^V P Al ) 

(^->|^5>) e -«<" } -<-' } )/ ft (20) 
{^)J MA f ^'(^<>|^>)e^"< 4> -^ } )/' i 

' l J -OO 

\ f ^ffi(0]*(<«M Na} >e" a ' ( ^' } " £ra,)/ *. 

' l J -oo 

where P(N a ) = ex.p[-/3(N a + l/2)hw a )/Z a is the prob- 
ability that there are N a phonons in the a th mode, and 
Z a = n ex -P(—P(N a + l/2)huj a ). The second and 
third lines come from current operator, cf . the second line 
of table [2] The fourth line of Eq.d20l) is the transition am- 
plitude induced by the transfer integral J, cf. the sixth 
line of table The fifth line of Eq.(l20b is the cc of the 
transition amplitude induced by the external field, cf. the 
second line of table Because we transformed the vi- 
brational displacements {xj} to normal coordinates {6 a }, 

the integral over displacements /[JJ dxj] in Eq.([T7l) be- 

j=i 

3Af 

comes an integral over normal coordinates /[ JJ dO a ] in 

(&^' a} \&{^' } ). Because A 1 is the final state of the LL 
transition induced by the transfer integral J, As is the fi- 
nal state of the LL transition induced by the external field, 
states A\ and As must not be occupied: one has factors 
[1 — f(EA 3 )] and [1 — /(Ea 3 )]- According to perturbation 
theory, one should (1) sum over all intermediate states; (2) 
sum over all the components of the final state. Because we 
don't know what the initial state is, we average over all 
possible initial states. At t = — oo, the external field and 
various residual interactions ( J, J, K and K') are not yet 
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Table 3 Transition amplitudes induced by the residual interactions and external field 



o y inuui 




A field A9 


ft J-00 L A 2 A V L /J \^A 2 1 A / c 


A field Bi 


ft J — 00 l B 1 A \ L )\ \^B 1 \*A / c 


B field A 9 




B field Bi 

— * 


4 /« ^q^(^o]*(4 <} i4^ } )e- lt,, '<" } - 4iVa})/R 

Al J — OO L ±J 1 ±J \ / J \ -D 1 1 -D ' 


AJ^A 3 


-i^,.<<»" , i*r- > >/!„*'="' < <°" , - 4 ''° ,,/ " 


A 4B3 




B ^ A 3 


-i(<"' > i^ 3B iHr- i >/:^*'e"'<" , -^ ,)/ » 


B ^B 3 





turned on, and the system is in equilibrium. The probabil- 
ities of \A • • • N a • • • ) is f(E A ) Y[ P(N a ). The first sum 

a 

comes from averaging over various initial states, the sec- 
ond sum comes from summing over the component of final 
state resulted from transfer integral, the third sum comes 
from summing over the component of final state caused by 
external field. The first factor in come from Eq.JTTl). One 
can similarly find the expressions of macroscopic density 
for other conduction processes. 

In an electric field E = Eo cos cut, the real part of 
conductivity tensor Rea a p can be extracted from the com- 
ponent of j a = (JapEp with time factor coscjt, the 
imaginary part Ima^ is extracted from j a with time fac- 
tor sin out. For the conduction processes in which the ini- 
tial state is a localized state, the expressions of Re cr a /3(uj) 
and Im a a p (qj) (a, f3 = x,y,z) are given in tabled the 
real part takes the upper sign, the imaginary part takes the 
lower sign. We used 

W A 3 A 2 = J drdp M (v - R Aa )^-0* 2 (r - R Aa ), 

(21) 

and 

v a 2 a 3 = j dr^ 2 (r - R A2 )^-0A 3 ( r ~ Ra 3 ), 

(22) 

etc to denote the matrix elements of velocity operator be- 
tween single particle states. One can show that the integral 

SJV 

over normal coordinates /[ JJ d6 a ], i.e. the sum over fi- 

a=l 

nal phonon states and averaging over initial phonon state 
can be carried out for all the contributions to the current 
densityj. The vibrational part of the current density j for 



each process becomes a time integral, cf. Sec. 2 of Ref. l22l . 
The time integrals I s , Q s and S s are functions of tempera- 
ture T and the frequency u of the external field, are listed 
in in Sec. 3.1 of Ref.|22|- Recr a p(uj) and lm(j a p((ju) sat- 
isfy the Kramers-Kronig relation, The reason that we list 
both of them is to show that (1) they are symmetric about 
the positive and negative frequency time integrals and (2) 
lm^(0) = 0. 

The conductivity from LL transitions fT3l derives from 
the processes depicted in Fig.[2al[3al|4a|and|4bl For small 
polarons, the Kubo formula was used in both Ref.Q and 
Ref.171 to obtain the second and fourth terms in table ffl 
but not the eighth and ninth terms. The reason is the fol- 
lowing. If one views both J and jf zeld as small param- 
eters, the eighth and ninth terms in table |4] result from a 
second order change in state, one in J and one in 
To obtain conductivity in Kubo formulation, the change in 
the density matrix of system is computed to first order in 
external field. By substituting this first order change into 
the macroscopic current density, one obtains the conduc- 
tivity by factoring out the external field |[T1. In this way, the 
new combinations of J and jf zeld (Fig|4a|and|4b]) are ex- 
cluded. In the present work, we apply linear response to 
external field at the last step, so that there are various com- 
binations between S and jf ield . From Eqs.(21, 25,27) of 
Ref. l22l . we see that the (cj, T) dependence of the fourth 
term is different from those of the eighth and ninth terms 
in table IH It will be interesting to see if new features of the 
eighth and ninth terms can be observed experimentally. 

If the initial state is an extended state £ B , substituting 
the corresponding ip^°\t) and ^ x \t) into the first term 
of Eq.([T7b, one can similarly derive the corresponding re- 
lation between the contribution to j and the diagram of a 
conduction process, see tables [2] and [3 There are 14 pro- 
cesses contributing to the conductivity, cf. Fig(6j(9} How- 
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Table 4 The conduction processes and conductivity: initial state is a localized state <p A . 



diagram conductivity 



|2^| 2Q S ^AA 2 ±111 y UJ AA 2 U A 2 AJ«\^A J ~'A 2 ) \ V A 2 AJ 
[I A ,A + ±lA 2 A-][l-f(E A Jf(E A ) 

m -^T,ABj^A Bl -^A)i{E\-ElX\<Ar 

[I B , A+ ±I Bl A-][l-f(E B ,)]f(E A ) 

133 -^^^A^A^A^AsAM^A-EXrHv^rjA.A 

[I A , A , A+ ±I A , AlA -]f(E° A )[l - f(E° Al )][! - /«)] 

-JKkT.AA 3Bl Im(w a A3Bl -v' Bl A 3 )(E A -E% 1 )- 1 (v^ A rj A3A 
^MfMglMJ! ~ Mjl W ~ K E A,)]f( E A) 

T.AA 2 B 3 Im ( W B 3 A 2 - V A 2 B 3 )(E A -Ea)' 1 (V A A )* J' B3A 



[33] 2hn s ^AA 2 B 3 \ UJ B 3 A 2 U A 2 B 3 )\^A 2 ^ A) \" A 2 j 

[I bm+ ±Ib 3 a 2 a-][1 - f(E B J[l - f(E A J]f(E A ) 



^T,ab 1 b 3 ^(< 3 b-^b 1 b 3 )J'b 3 a(K-Eb 1 ) 



_1 r^ v 

F^Jl 2hS2 s *-~iA±l 1 ±l 3 V — B 3 B X ~D\D 3 J~ B 3 A\~ A ~D\J \ B\A) 

[Ib^a + ^b 3Bi a-][1 ~ f(E B J][l - f(E B ,)]f(E A ) 

+^J2 A2AlA ^(w AA2 -v A2A )(E Ai -E A2 y 1 (v a A2Ai yj AlA 

(QlAzAj A+^ z QlA 2 A 1 A- 

m +lEk ^A 2 A lA I™ «A 2 -Va 2 a)Ja 2Ai {E A -E Al )-\v AiA y 

(Q2A,A,A + ±Q2A*A,A-)f(E A )ll ~ f( E Ao )] 

^ +^Y.B 2 A 1 A^(<B 2 -< 2 A)(E Al -E% 2 y\v% 2Ai yj AlA 

(Q 1b ,a,a + ±Qib*a,a-MEa)11 ~ f(E B J] 



+ §k^A 2 A 1 A^(<B 2 -^kA)JtAS E A-EXr\< 1 AT 

)f(E A )[l-f(E B )] 

+ %&ZA 2 B lA ^(<A 2 -VA 2 A)(E A -E B J-\v BiA y 
(Q?L B ,A + ±Q?MB,A-)f(E A )[l ~ f(E A J] 
+ §k^A 2 A 1 A^(<A 2 -< 2 A)(E% 1 -E%y\v A2Bi yj'£ lA 

)f(E A )[l-f(E A )] 

+ ^J2B 2 B 1 A^(<B 2 -< 2 A)(E A -E° Bl )- 1 (v a BiA y 
(Q? B ? B ,A + ±Q?B*B,A-)f(E A )[l ~ f(E B J] 
+ ^J2B 2 B 1 A^(<B 2 -^kA)(El 1 -E% 2 y\v% 2Bi yj' B * lA 

(Q2A 9 Bi A+^ Z Q2A 2 B 1 A- 

)f(E A )[l-f(E A J] 



ever three order K 1 EE transition processes (Fig |7dl[9cl[9dl) 
are zero due to the e-ph interaction selection rule. For pro- 
cesses involving only extended states, order K 2 processes 
are the first nonzero contributions. The first term in Eq.JTTt 
indicates there are 4 such terms (FigfTQl). Thus if the initial 
state is an extended state £ B , to lowest order self-consistent 
approximation, there are 15 conduction processes. The cor- 
responding expressions for conductivity are listed in table 
The time integrals are given in Sec. 3.2 of Ref. ll22l . 

For EE transition, the order K 1 contributions to con- 
duction are zero (Fig|7dl[9cJ[9d|), so that one has to take into 
account order K 2 contribution (Fig[lOa| [TObl [TOcl [TOdb . 
The time integrals are given in Sec. 3.3 of Ref. l22l . The 



conductivity from the pure EE transitions derives from the 
second (Fig(6bJ carriers are scattering by the static dis- 
order), and the thirteenth to the sixteenth terms (FigfTOl 
the carriers in extended states are inelastically scattered by 
phonons) in table [5j This is consistent with the Boltzmann 
theory: in the lowest order approximation, the scattering 
probability is second order in the e-ph coupling constant K 
while the distribution function is order K° (the non- 
interacting elementary excitations) 1 26]. Thus the collision 
integral is proportional to K 2 and the change in dis- 
tribution function and the resulting conductivity is propor- 
tional to K 2 . 
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Table 5 The conduction processes and conductivity: initial state is an extended state £ B . 



diagram 






conductivity 


M 






-^T, BA ^i{E\ 2 -E%)-\ BA2 { WBA2 -vX B ) 
\I A ,b + (") ±^,b-M][1 " f& A Jf(E B ) 


m -%£Y.b Bi lrai{E%-El)- l vi Bi {w% B -v% lB ) 

[I BlB+ M ±7 BiB _H][l - f(E B J]f(E B ) 


M 


2nn s L-i 


BA 3 A 2 lm ( E A 2 - E B)~ 1 V BA2 (w A3A2 -V A2A3 )[I A3A2B+ (L0) ± I ^^{u)} 

[l-/(^J][l-/(£? A ,)]/(£? fl ) 




N e e 2 
2hfi s 


E 


BBi^s 1111 ( £; b 1 -- E b)" 1 v BBi (^ 3Bi -v BiA3 )[/_ 43BiB+ (o;) ±7 A;jBiB _(w)] 
[1-/(£ A J][1-/(£ B J]/(£ B ) 


M 


Nee- 1 

2nn s 


E 


BB3A2 ImK 2 - J BB)- 1 ^ 2 (^B3A 2 -^ 2 B 3 )[^ 2 B3B + H±^ 2 B3B-H] 

[1-/(S A ,)][1-/(S B ,)]/(£? B ) 


M 






+ Ea 2 a,b Im K A2 - V S 2 b)« -^ 2 )" 1 «a 1 )* 

(^^B+i^A.B-)/^)^ - /(^A,)l 


[8bl 






+ I^EA 2 A 1 Blm(^ 2 -^ 2 B)^ 2 A 1 (^B-SA 1 )- 1 (< i B)* 
( 5 '2A 2 A 1 B+ ±5 '2A 2 A 1 B-)/(-E B )[l - 


[8cl 






+ MEa 2 a 1 bIhi( WbB2 -^b 2 b)^ 224i K-^°a 1 )- 1 (<b)* 

( 5, 1B,A 1 B+ ±5 'iB2AiB-)/(-Eb)[ 1 ~ f( E B,)] 


[8dl 






+iS^A lB imK B r^B)(^ 1 -4 2 )- 1 (< 2Al )' 

( 5, 2B,A 1 B+ ±5 '2B 2 A 1 B-)/(-Sb)[ 1 ~ f( E B,)] 


M 






+ ^T,B 2 A 1 B^^A 2 -vXB)( E % 1 - E A 2 y\v A2 B i r 

(5f A , BiB+ ±5f A2BlB _)/(£; B )[l - 


M 






+ ^E^ 2 B 1 Blm(^BA 2 -«A 2 B)(SB-^B 1 )- 1 (<B)* 
( 5, 2A 2 B 1 B+ ±5 '2A 2 B 1 B-)/(-Sb)[ 1 ~ /(^Jl 


fTOal 






+ 2^ Eb Bi B 2 B 3 toK 3 Br"B 1 B3)i(4-£ B 1 )- 1 
(^bH^B^B+H ± ^B^B-M] 

[l-f(E Bi )][l-f(E B Jf(E B ) 


[TObl 






-^E B3B2BlB imK B3 -^ B )i(4-4 1 )- 1 « B )* 

(• ft 'lB 2 B,B 1 B+ ± - ft 'lB3B 2 B l B-)/(-E ; B)[ 1 _ /C^B,)] 


HOcI 






( K 2B^B 9 Bi B+ ±K 2B 3 B 2 B 1 B-)f(E B )[l - f{E B J\ 


flOdl 






-^E^^s^Im^-^sW^-^)" 1 ^^)* 

(- FS: 3B 2 B,B 1 B+ ± ^3B3B 2 B 1 B-)/(£'b)[ 1 ~ /("^bJ] 



(°)W,( 2 )* 



The expression for the conductivity is a sum of 29 
terms listed in tables S] and [3 



Table 6 Contributions to a yx from V^ (1) * 



J° 




Figl 


lllbj 




J 




Figl 


[Hal Fid 


[T2bl 


J 




Figl 


ESI Fig| 


H2dl 



Table 7 Contribution to from 

^( 2 )*v^ (0) ) 



J° 


< v^r 




FigfTTal 


J 






FigfTIal 


J 






Fisll2blll2clll2dl 



6 Hall mobility If an amorphous semiconductor is 
placed in an external magnetic field, the conductivity is 
still defined through the current density j a = J2b a a/3E/3- 
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M.-L. Zhang and D. A. Drabold: 




a) b) 

Figure 6 Initial state is extended state order [K']° and 
order K° contributions to conductivity 




Figure 7 Initial state is extended state £ B : order [K'} 1 and 
order K 1 contributions to conductivity 




a) b) c) d) 

Figure 8 Initial state is extended state £ B : orderiT, [J'} 1 
and order J 1 contributions to conductivity 



If we apply a static electric field along the x direction and a 
static magnetic field along the z direction, the Hall voltage 
is along the y direction and proportional to both B z and 
E x . Amorphous semiconductors are isotropic, and the Hall 
mobility is given byl8l: 

= B~ 1 a yx /(j xx . (23) 




a) b) c) d) 



Figure 9 Initial state is extended state \ order [K'\ l and 
order K 1 contributions to conductivity 




Figure 10 Initial state is extended state \ order K 2 con- 
tribution to conductivity 

The non-diagonal conductivity in a magnetic field is de- 
scribed by the second term of Eq.(fT7b. 

In general, both the carriers in localized states and 
carriers in extended states contribute to the Hall effect. 
There are too many terms in the full expression for non- 
diagonal conductivity to reproduce here. We restrict our- 
selves to the LL transitions. This is a reasonable approxi- 
mation for intrinsic and lightly-doped amorphous semicon- 
ductors, where carriers in extended states are rare. At the 
end of this Section, we will indicate how to obtain the full 
expression of Hall mobility. 

By substituting the corresponding ijj^ and 
into the second term of Eq.(fT7b, one finds the current 
density in a magnetic field and an electric field. Various 
contributions may be visualized by diagrams as for the 
conductivity: the transition amplitudes in Table [3] are still 
applicable. The expression for the two components con- 
nected by the current operator (^ (0) V^ (2) * - ^ (2) * V^ (0) ) 
is the same as that for two components connected by 
(V^ ^ V?/^ 1 ^* — V?/^ ^) in the ordinary conductiv- 
ity. We still use a wavy line to depict such a expression: 
it points from one component of ip^ to one component 
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of ^ (2) *. There is one order J u process and four order 
J 1 processes. In Table |7] we list the contributions from 
(^ (0) V^ (2) * - ^ (2) *V^ (0) ). If we exchange ip {0) and 
?// 2 )*, the current density changes sign. The new contribu- 
tion to current density from t/j^ V?/^ 1 ^* does not include 
this exchange antisymmetry. We represent it by a dashed 
arrow line, cf . the last line of table [2j The dashed arrow 
line points from one component in ijP^ to one component 
m There is one order J u process and four order J 1 

processes. We list the contributions from V?/^ 1 ^* in 
Table [6j Order J° contributions are illustrated in FigfTTT 
Order J contributions are illustrated in Fig[l2land Fig[T3l 
In the symmetric gauge A = |B x r, the coupling be- 
tween two localized states by the external magnetic field 
becomes 



J 



field _ 
A 2 A 1 - 



2m 



(24) 



where x M Al = J dr<f> A2 x4> Ai ,andLf 2Al = J dr4> A2 L z <p^ 
is the matrix element of the z component of electronic or- 
bital angular momentum. To compute the Hall mobility, 
we require only the terms which are proportional to B Z E X . 
For example in (jfcf 2 )* {J^f A \ )*, we only keep 

{e 2 /2m)B z E x {L z MA2 x MAl + x A3M L z MAi )* . (25) 




AC i A i 

l IT I 

a) b) c) d) 

Figure 12 Order J 1 contributions to a yx : I12al results 
from V'jo ) V{V' (1) [A 4 A 1 A^)*.fT2bl results from 
tp { j} [A J '-T Ai 4 A 2 ]}*.[l2c]results from ip {1) [A 
J -X Ai J/ 4' d A 2 ]V^^ [Ed] results from ^ (l) [A. J '^° 
^A 2 ]V^V*. 




b) 




d) 



Af j-jiBldj* 

U 




a) 



Figure 1 1 Order J° contributions to q^. lllal results from 



( )x7,/,( 2 )* 
J 



^jo*Wjo )-EHresults from ip^JVip 



,,(o) 



One can extract a yx from j y . The corresponding ex- 
pressions for various processes are given in table [8] The 
time integrals U, U', Z, Z' and Y are given in Sec. 5 of 
Ref. ll22l . they are functions of temperature T. The non- 



diagonal conductivity a yx is a sum of 10 terms listed in 
table El In the previous theory of the Hall effectEHU, 
researchers introduced magnetic field dependent Wannier 
functions as the basis and expanded the magnetic field de- 
pendent phase in the final step. Such a method excludes 
the interference between the transition amplitudes caused 
by electric field and magnetic field. In the present work, we 
treat electric field and magnetic field on the same footing. 



Figure 13 Order J contributions to <T,,^: ll3al results from 
VV'jo* - ^| 7 2 o ) *VV'^ ) ).[l3B[l3i]and[l3d|result from 



A mixing similar to Eq.([25t appears in each term of table 
El In the application of Kubo formula OEE, FigOIEEl 
I13b[ I13c[ I13dl are excluded. These new terms are of the 
same order as Fig |llal and ll3al but the temperature depen- 
dence are different. To check the existence of these new 
term, it requires more experimental data points lfT3l . 

The full Hall effect from the carriers in localized states 
and in extended states can be found by eventually chang- 
ing the states from A to B, and correspondingly changing 
the residual interaction and the transition amplitudes. For 
example there are 7 more diagrams for each diagram in 
FigEl Using tables [2] and one can easily write out the 
expressions for each diagram |[T3l . 

7 Conclusion The microscopic response method |[T3l 
OIH is applied to systematically estimate the conductivity 
(tables 4 and 5) and Hall coefficient (table 8 and Eq.(l23l)) 
in amorphous semiconductors. A diagrammatic represen- 
tation of transport coefficients is introduced by inspecting 
the structure of the observed current density and the per- 
turbation expansion of wave function about small param- 
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M.-L. Zhang and D. A. Drabold: 



Table 8 diagrams and ay 



diagram 



a 



yx 



TTal 



B ^2^0TF S T,aa 2Ai lm{UAA 2 Aj(E A )[l - f(E A 

l{Wy 



^-v^ A )(L^x AlA +x A2Al L^ A y} 



TTbl 



U'aa 1 a 2 (T)f(E A )[l - f(E Al )][l - f{E A2 )]} 



A 2 A\* 



12al 



B 



2mh 3 n s ^AA 1 A 2 A 3 
A 3 A 2 ( t Ax A 

y 



lmf(E A )[l-f(E A2 )][l-f(E Aa 



" 2 (L? 3A x* A2Al JX lA + x AaA (Lf 2Al )* J AlA )Z AAsA2Al 



12b1 



^2^fe EAAxA^a Im/(^A)[1 - f(E A2 W - f(EA 3 )] 

w A3A Hx A3A r A2Al (L AlA r + L A * A J* A2Al x* AlA )Z AA3A2Al 



12cl 



3^ 



1 + (^z 3 )*^A 2 A!) Ja 1 aZ / AAsA2Ai 



1251 



-s 2 



2mh 3 f2. 



E 



Im/(^)[l-/(^ 2 )][l-/(^ 3 )] 

{•^\ 3 A^A 2 AiL AlA + (^ 3A )Va 2 A 1 XA 1 a)^aA 3 A 2 A 1 



AA\A 2 A 3 



1551 



iV e e J 



E 



AA\A 2 A 3 



lmf(E A )[l-f(E A2 )}[l-f(E Aa )} 



2mh 3 Q s 

A 3 A 2 _ v f2A 3 )j AaA ( L A 2 A lxAiA + ^ A2Ai L^)*Z AA3A2Al 



he 



TT^ri^^z ^2a 3 A 2 A 1 A ^ m ^ r A 3 A 2 A 1 Af(EA)[i- ~ f(EA 3 )] 
iJl 3 A 2 (w V AA, ~ v A^A)[(-^A 2 A 1 T( x A 1 Ay + 



I3cl 



7 ^ fe B,^EA 3 A 2 A 1 Alni^3A 2 A 1 A/(£;A)[l-/(£;A3)] 



13dl 



fl B 2m¥^ EAsA.AiA^^sA^A/^A)!!-/^)] 
)[(^A 3 A 2 )*(^A 2 A 1 )* + (^A 3 A 2 )*(^ 2Ai ; 



^IxA^AAs "^A 3 A 



eters and about external field. From the topology of the 
diagrams, one can easily write down various contributions 
for a transport coefficient to a given order of small parame- 
ters. This is helpful for higher order processes. Comparing 
to the Kubo formula, all important conduction processes 
can be easily taken into account by the diagrammatic rules 
in tables [2] amd [3j Impressive advances have been made 
for computing transport properties from first principles us- 
ing the Keldysh formalism! 27 , 28 , 29 ] . The present method 
is complementary to that work as a means of computing 
temperature dependence of the transport coefficients rather 
than I-V curves. 

The present work has two obvious advantages over the 
thermal average method of classical MD: (i) the vibrations 
of atoms are described by quantum mechanics, the results 
obtained are valid for arbitrary temperature (at least for 
which the harmonic approximation is valid), whereas the 
results of MD are correct only above the Debye temper- 
ature of a material; (ii) In the expressions for conductiv- 
ity and Hall mobility, the e-ph interaction is already inte- 
grated out. To calculate the transport coefficients at a given 
temperature, one only needs the eigenvectors and eigen- 
values of the Kohn-Sham Hamiltonian and the dynamical 



matrix for one configuration. One does not need to com- 
pute the time-consuming MD trajectory or worry about 
proper equilibration. The method is ideally suited for ab 
initio simulation and is in the process of numerical imple- 
mentation. 
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References 

[1] R. Kubo, J. Phys. Soc. Jpn. 12, 570 (1957). 

[2] D. A. Greenwood, Proc. Phys. Soc. (London) 71, 585 (1958). 

[3] M.-L. Zhang and D. A. Drabold, Phys Rev. E.83, 012103 

(2011). 

[4] M.-L. Zhang and D. A. Drabold, fcarXiv: 101 1.4 11 3J/1 [cond- 

mat.stat-mech], J. Phys: Condensed Matter, in press. 
[5] I. G. Lang and Yu. A. Firsov, Zh. Eksperim i Teor. Fiz. 43, 

1843 (1962) [Soviet Phys.-JETP 16, 1301 (1963)]. 
[6] Yu. A. Firsov, Fiz. Tverd. Tela 5, 2149 (1963) [Soviet Phys.- 

Solid State 5, 1566 (1964)]. 
[7] J. Schnakenberg, Z. Physik. 185, 123 (1965). 
[8] T. Holstein and L. Friedman, Phys. Rev. 165, 1019 (1968). 
[9] A. Aldea, J. Phys. C: Solid State Phys. 5, 1649 (1972). 



Copyright line will be provided by the publisher 



pss header will be provided by the publisher 



13 



[10] S. K. Lyo, J. Phys. C: Solid State Phys. 6, 2158 (1973). 

[11] N. F. Mott and E. A. Davis, Electronic Processes in Non- 
crystalline Materials, Clarendon Press, Oxford, (1971). 

[12] M.-L. Zhang and D. A. Drabold, Eur. Phys. J. B. 77, 7-23, 
(2010). 

[13] M.-L. Zhang and D. A. Drabold, Phys. Rev. Lett. 105, 

186602 (2010). 
[14] R. A. Marcus, Rev. Mod. Phys. 65, 599610 (1993). 
[15] T. Holstein, Ann. Phys. 8, 325 (1959); 343 (1959). 
[16] D. Emin, Phys. Rev. Lett. 32, 303 (1974); Adv. Phys. 24, 305 

(1975); Phys. Rev. B 43, 11720 (1991). 
[17] P. A. Fedders, D. A. Drabold and S. Nakhmanson, Phys. Rev. 

B58, 15624 (1998). 
[18] A. Miller and E. Abrahams, Phys. Rev. 120, 745 (1960). 
[19] P. Thomas, J. Non-Cryst. Sol., 77&78, 121 (1985). 
[20] D. A. Drabold, P. A. Fedders, S. Klemm and O. F. Sankey, 

Phys. Rev. Lett, 67, 2179 (1991). 
[21] R. Atta-Fynn, P. Biswas and D. A. Drabold, Phys. Rev. B69, 

245204 (2004). 

[22] see supporting information at 

http ://onlinelibrary. wiley. com/store/ 

[23] C. M. Van Vliet, Equilibrium and Non-equilibrium Statisti- 
cal Mechanics, World Scientific, Singapore (2008). 

[24] M.-L. Zhang and D. A. Drabold, Phys. Rev. B81, 085210 
(2010). 

[25] D. Bohm, Quantum Theory, Prentice-Hall, New Jersey 
(1951). 

[26] E.M. Lifshitz and L. P. Pitaevskii, Physical Kinetics, 

Butterworth-Heinemann, Oxford (1981). 
[27] M. Brandbyge, N. Kobayashi, and M. Tsukada, Phys. Rev. 

B 60, 17064 (1999). 
[28] M. Brandbyge, K. Stokbro, J. Taylor, J. L. Mozos, P. Orde- 

jon, MRS Conf. Proc. 636, D9.25 (2000). 
[29] J. Taylor, J. Wang and H. Guo, Phys. Rev. B, 63, 245407 

(2001). 

[30] D. Binosi and L. TheuBl, Computer Physics Communica- 
tions, 161, 76 (2004). 



Copyright line will be provided by the publisher 



